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HEAT KERNEL ESTIMATES FOR THE BESSEL DIFFERENTIAL 

OPERATOR IN HALF-LINE 
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Abstract. In the paper we consider the Bessel differential operator + 

dx‘‘ X ax 

in half-line (a,oo), a > 0, and its Dirichlet heat kernel {t,x,y). For ^ = 0, by combining 
analytical and probabilistic methods, we provide sharp two-sided estimates of the heat kernel 
for the whole range of the space parameters x,y > a and every t > 0, which complements the 
recent results given in [1], where the case y ^ 0 was considered. 


1. Introduction 

We consider the Bessel differential operator 

T id ) = I 2 /x -h 1 d 


// G R, 


X dx ’ 

on the half-lines (a, oo), for every a > 0. Let us denote by p ^\ t , x , y ) the corresponding Dirichlet 
heat kernel considered with respect to the reference measure m^^\dy) = For y = 0 

we will simply write Pa { t , x , y ) and we will also omit superscripts in the notation of the other 
objects related to that case. The explicit formula for p ^^\ t , x , y ) := p ^\ t , x , y ) in terms of the 
modified Bessel functions is well-known 


P^^Ht,x,y) = ^{xy) ^exp(- 


1 


-h y2 

Yt 




t ) 


x , y , t > 0 . 


For a > 0 and y 7 ^ 0 the sharp two-sided estimates of pi^^ { t , x , y ) of the form 


pi^Ht,x,y) 


1 A 


(x - a)(y - a) 


t 


1 A 


xy\\d\-^ 1 1 

t ) (xy)/^+^/2 ^ 


exp 


(x - yf 

2t 


( 1 . 1 ) 


( 1 . 2 ) 


r^n 

for every x,y > a and t > 0, were obtained recently in [T]. Here ~ means that the ratio of 
the functions on the left-hand side and the right-hand side is bounded from below and above 
by positive constants (depending only on y) for the indicated range of variables. Notice that 
the exponential behavior is here described explicitly, i.e. there are no different constants in the 
exponential terms in the lower and upper bounds. This is also the case with the main result of 
the present paper, where the two-sided sharp estimates of pi{t,x,y) (for y = 0 ) for the whole 
range of the space and time parameters are provided in 


Theorem 1. ITe have 

Pi { t , x , y ) « In X In y ( In 

and 


3t 


: In ■ 


3t 


pi{t,x,y) ^ A 
for every x, y > 1 and t > 0. 


X + y + Vi 

(x-l)(y-l)\ 1 


-exp 


x^ + y^ 
2t 


xy < t. 


t 


VYJi 


exp 


{x - yf 

2t 


xy > t, 


(1.3) 


(1.4) 


2010 Mathematics Subject Classification. 35K08, 60J60. 

Key words and phrases. Bessel differential operator, heat kernel, sharp estimate, half-line, Bessel process, 
transition probability density. 

The project was funded by the National Science Centre grant no. 2013/11/D/ST1/02622. 

1 

















2 


KAMIL BOGUS, JACEK MALECKI 


Using the scaling property 

Pa{t,x,y) = -pi{t/a^,x/a,y/a), x,y > a, t > 0, 
a 

one can immediately deduce the corresponding result for any a > 0. Moreover, the result can 
be rewritten in one formula covering both cases, i.e. xy/t small and xy/t large. More precisely, 
taking into account the asymptotic behavior of Io{z) at zero and infinity (see (I2.2p l. we can 
write the estimate in the following way 


Pa{t,x,y) 

P{t,x,y) 


1 A 


ln(a:/a) ln(y/a) 


3xy + 3t 3xy + 3t 

In- -= In- jL 

ax + ay t ay + ov t 


-1 


xy 

i + T 


for every x,y > a and t > 0. Note that the constants appearing in the estimate are absolute and 
do not even depend on a. Notice that for xy > t (and a = 1) the estimates in (|1.2H are exactly 
of the same form as those given in m- However, for xy < t the logarithmic terms appear in 
which is not the case in m- Moreover, the logarithmic behaviour of Pa{t,x,y) require 
more delicate methods, i.e. those used in [1] can not be applied to prove m- 

Since exploring the behaviour of the heat kernels is very important in many areas of Math¬ 
ematics (for example in differential equations, harmonic analysis, potential theory or theory of 
stochastic processes), the literature abounds with a huge amount of different kinds of deep and 
beautiful results on heat kernels estimates, which are difficult to briefly describe. However, even 
if we deal with the most classical case, i.e. the heat equation based on Laplacian in R” (see 
m or the works of E.B. Davies m, m, 0,0) or we consider the corresponding problem on 
manifolds (see for example HH and the references within) or finally, we look at the resent results 
on the Fourier-Bessel heat kernel estimates on (0,1) studied in |12| and |13) . the overwhelming 
majority of the results are qualitatively sharp, i.e. the constants appearing in the exponential 
terms in the upper and lower bounds are different. The presented result given in Theorem [H as 
well as the result of [T], is free from this defect and here the exponential behaviour is described 
explicitly. However, this requirement imposes more complicated form of the part describing 
non-exponential behaviour of p{t,x,y). We refer Reader to the Introduction in [T], where more 
detailed discussion on this phenomenon was conducted. 

The paper is organized as follows. The basic properties of the Bessel function Io{z) as well 
as the notation related to Bessel processes are recalled in Preliminaries. Moreover, in Lemma [H 
we introduce some technical estimates of a derivative of the heat kernel on (0, oo), which will be 
useful in the sequel. Then, in Section |3l we provide the proof of formula (11.311 which is divided 
into two parts. The first one relates to the case, when x‘^/t and y‘^/t are bounded (Proposition [1]) 
and in Proposition [2] we deal with the case when (x V y)‘^/t is large. Finally, in Section 0] we 
justify ([HD. 


2. Preliminaries 


2.1. Modified Bessel functions. We begin with collecting some selected properties of the 
modified Bessel function Iq{z) used in the sequel. We start with the following definition 


Io{z) = 

71=0 




z eK. 


The above given series representation immediately implies that Io{z) is an increasing function 
on [0, oo). On the other side, the function e^I(){z) is decreasing on [0, oo), i.e. we have (see Thm. 
2.1 in (H) 

^ y>x>0. (2.1) 

Io{x) 

Obviously /o(0) = 1 and the behavior at infinity is described by (see |8j 7.13.1 (5)) 

7o(z)~^(l + 0(l/z)), 
y27r2; 


z ^ oo. 


( 2 . 2 ) 
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Finally, we recall the estimates of the ratio of Io{z) and its derivative Ii{z) given in |11| 

z It{z) 2z 


h{z:) 

< < 


z + 2 Io{z) 2z + 1’’ 


z > 0. 


(2.3) 


2.2. Bessel processes. We write and for the probability law and the corresponding 
expected value of the Bessel process with index /i > 0 on the canonical path space with 
starting point = x > 0. The filtration of the coordinate process is denoted by = 

a{R^t'^ : s < t}. The transition probability density (with respect to the reference measure 
m{dy) = ydy) of two dimensional Bessel process {y. = 0) is given by 


P{t,x,y) = ^exp ( - 


+ 2 /^ 

¥t 


^ (y) , x,y,t>0. 


(2.4) 


The laws of Bessel processes with different indices are absolutely continuous. In particular, the 
corresponding Radon-Nikodym derivative for the laws of the Bessel processes with index 0 and 
1/2 is described by 


dP^°'> 


dP 


( 1 / 2 ) 




( 1 / 2 ) 


wit) 


- 1/2 


exp I I 


ds 


8 Jo w‘^is)J ’ 


(2.5) 


where x > 0, and the above given formula holds Pi^'^^^-a.s.. 

Continuity of the paths and the Hunt formula enable us to write the corresponding transition 
density function for the process r[^^ killed when it leaves a half-line (o, oo), a > 0, in the 
following way 


Pait,x,y) = pit,x,y) - qx,ais)pit - s,l,y)ds, x,y > a, t > 0. 

Jo 


( 2 . 6 ) 


Here qx,a{s) denotes the density function of Tq, the first hitting time of a level a by the process 
starting from x > a. Its uniform sharp estimates were provided in [3] (see Theorem 2 therein). 
More precisely, we have 


fe,i(s) 


X — 1 1 


3/2 


and 


Qx,i{s) 


X — 1 


X s 


1 -|- Inx 


exp 


(x — 1)" 

2 ^ 


1 

-- exp 


s < 2x, 


(x- 1)^ 
2s 


s > 2x. 


(2.7) 


( 2 . 8 ) 


X (1 -I- ln(s-|-x))(l -I- ln(l-|-s/x)) s 
These estimates lead to the uniform bounds of the survival probabilities (see Theorem 10 in m) 

Inx 


qx,iis) ds « 1 A 


ln(l -I- Vt) ’ 


X > 1, t > 0. 


Note also that 


/ OO 1*00 

pi{t,x,y)m{dy) = BxiTi > t] = J qx,iis)ds. 


(2.9) 


( 2 . 10 ) 


We end this section with the following very useful lemma, which refers to the uniform estimates 
of the derivative in x of pit, x, y). 


Lemma 1. ITe have 


dx 


pit,x,y) pzpit,x,y) ■ X 


whenever xy <t and > 4f. 
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Proof. Differentiating the right-hand side of (|2.4I) and using Iq{z) = Ii{z) lead to 

for every x, y > 0 and t > 0. The upper bounds in (I2.3p gives 

h{xy/t) 2xy 2xy‘^ 

y-rz - rr — x <y ■ - x < -. 

Io{xy/t) 2xy + t t 

On the other-side, using the lower bounds in (12.311 and the conditions xy <t,y‘^> At, we get 

h{xy/t) xy^ ^ xy"^ 1 xy^ 

^Io{xy/t) ^-xy + 2t ^ “ 3t ^“12’ t' 

This ends the proof. □ 

Remark 1. If we assume that 2x < y and y^ > 4t (namely, we replace the condition xy < t by 
2x < y) we can write, using lower bound given in (12.311 . that 

h{xy/t) xy2 x{y‘^-xy-2t) x{y‘^/2-2t) 

y-—- tt — X > -X =- >- > 0 

Io{xy/t) xy + 2t xy + 2t xy + 2t 

and, as a consequence, we obtain that the function x i—)• p{t, x, y) is increasing for y^ > At and 
1 < X < y/2. 


3. Estimates for xy/t small 


By symmetry, it is enough to consider 1 < x < y. Moreover, in this section we only deal with 
the case xy < t. ft implies that x'^ < xy < t and consequently m reads as 

In X In y , 

pi{t,x,y) ^—^—p{t,x,y), y <t, 

In 3t 

and 

, . In X , . 2 

pi[t,x,y) Tn j-^^^p{t,x,y), y >t. 

ft naturally splits the proof into to parts related to the case when y^/t is small (Proposition [1]) 
and when y^/t is large (Proposition [2]). 


Proposition 1. For every m >2 we have 


m Inxlny 

pi{t,x,y) ^ ^ p{t,x,y), 

In (3t) 


whenever xy < t and y^ < mt. 


Proof. Taking x ^ 0"*“ in (j2.1ll we get that loiz) < for every z > 0, which together with the 
obvious general estimate pi{t,x,y) < p{t,x,y), see ( 12 . 6 p . gives 


Piit,x,y) < - exp - 


t 


x^ -|- y^ 

Yt 


Io{xy/t)<-, x,y>l,t>0. 


Applying this estimate to the middle term in the Chapman-Kolmogorov equation given below 
we arrive at 

/ oo poo 

J pi{t/3, X, w)pi {t/3, w, z)pi {t/3, z, y)m{dw)m{dz) 

2 POO POO 

<-y pi{t/3,x,w) m{dw) J pi{t/3,y,z)m{dz) 


m Inx Iny 


0 POO POO 


where the last line follows from the fact that p{t, x, y) ~ 1/t, whenever xy < y'^ < mt, and (|2.9p 
(note that t > xy > 1 and consequently ln(l + y^t/3) ~ ln3t). 

















HEAT KERNEL ESTIMATES FOR THE BESSEL DIFFERENTIAL OPERATOR IN HALF-LINE 


5 


To deal with the lower bounds we begin with showing that there exists a constant Cm > 0 
such that 


pi{t,x,w) > c„ 


Inx 

ln(3t) 


pit,x,w) 


(3.1) 


whenever 1 < x < w(2 and re > 1 + 2^/t. Indeed, denoting the subtrahend in (12. 6 p as ri(t, x, y), 
we can write 


Pi {t, x,w) _ ^ _ n {t, X, w) p{t, 1, w) 


pit,x,w) 


where using (EH) we have 

i^,x,w) 

p{t 


.,x,w) f 
-j-^ = / 9x,i(s 

A,w) Jo 

/ 9x,i(s) 
Jo 


t 


t — s 
t 

t — s 


exp 


exp 


p{t,l,w) pit,x,w) 

w'^ + lf 1 l\\ lo {w/{t - s)) 

{w - 1)2 


t — S t J J lQ{w/t) 
1 1 ' 


ds 


t — s t 


ds. 


Consequently, using monotonicity of the function u i-A u~^ exp [—{w — l)‘^/2u) for tc > 1 + 2\/t 
and 0 < u < t, we get 

rt 


ri{t,x,w) 


< / qx,iis)ds. 


p{tA,w) Jo 

On the other hand, the monotonicity of x —>■ p{t,x,w) showed in Remark [1] gives p{t,l,w) < 
p{t,x,w). Hence, using (12.9p . we obtain 


pi{t,x,w) 


pt poo 

>1- / qx,ii^)ds = / qx,i{s)ds 

Jo Jt 


p{t,x,w) 

Finally, the Chapman-Kolmogorov equation and (13.ip give 


Inx 
In 3t 


^ Cm, 


Since Io{z) > /o(0) = 1; we get 


pZy/mt 

Pi{t,x,y)> / pi{t/2,x,w)pi{t/2,y,w)m{dw) 

J 2\/mt 

In X In y 


ln2(3t) J2./^t 


p(t/2, X, w)p{t/2, y, w) m{dw) 


Inxlny 1 

Pi{t,x,y) > Cm 


x 2 y 2 


pSy/nii 
12\/mt 


exp 


2w^ 


m{dw) 


m Inxlny 1 m Inxlny 

m{dw) ss p{t,x,y), 


ln2(3t) t2 J^^^ 


ln2(3t) 


where the last estimate once again follows from p{t,x,y) ~ 1/t, which is true whenever xy < 
y2 < rnt. The proof is complete. □ 


Proposition 2. ITe have 


Pi{t,x,y) 


Inx 


■Pit,x,y), 


ln{3t/y) 

whenever xy < t and > \Qt. 

Proof. Note that for y < 52 the conclusion follows from Proposition [1] with m = 52^. Thus from 
now we will additionally assume that y > 52. The second simplification comes from the fact 
that p{t,x,y) ~ p{t,l,y) whenever xy < t. Taking that into consideration and using (|2.6I) . we 
can write 

pi{t,x,y) p{t,x,y)-p{t,l,y) , x , , p{t,l,y) - p{t - s,l,y) 


p{t,x,y) 


p{t,i,y) 


poo pt 

+ / qx,i{s)ds + 

Jt Jo 


p{tA,y) 


-qx,i{s) ds. 
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To deal with the first part, we use the Mean-value theorem and Lemma [T] to show that there 
exists 1 < Xc < X such that for xy < t and > 16t > At we have 

,,2 


Q 

P{t, X, y) - p{t, 1, y) = (x - 1) • X, y) 


Xc{x - l)^p(t,Xc,y). 


In particular, it implies that the difference is positive. Moreover, since p(t,Xc,y) ~ p{t,l,y), 
there exists constant ci > 0 such that 

p{t,x,y) - pit, l,y) < Cl • x(x - 1) ^p{t,l,y). (3.2) 

The estimates of the second part are given in (12.9p . Finally, we split the last integral into 
three parts 



p(L ^,y) -pjt - 
pitA,y) 


s, i,y) 


qx,i{s)ds 


^2 -h ^3 


and estimate each of them separately by examine the behavior of the ration appearing under 
the integrals for the corresponding range of parameter of integration. More precisely, using the 
inequality e~^ >1 — 2 and the fact that loiz) ~ 1 near zero we get the existence of the positive 
constant C 2 such that for every s < At'^/y'^ < t/A we have 

P(L1,1/) -Pit-s,l,y) ^_ t_ ^oiv/jt - g)) 

p{t, 1, y) t-s V -s) J loiy/t) 

< ^_ t_ A _ (1 + loiy/jt - s)) 

t- s\ 2t{t - s) J hiy/t) 

_ t Ipiy/jt-s)) jl + y‘^)s lojy / jt - s)) y^ 

t-s hiy/t) 2 (f-s )2 hiy/t) -'^ 2 ^ 2 ^- 

Note that the difference of the first two components is negative and has been omitted in the 
last estimate. Since s < At'^/y'^ < t/A, we have 2 t{t-s) — ~ ^/‘A- Thus, the inequality 

e~^ < 1 — 2 / 3 , which is true for 0 < 2 < 8/3, gives 

Pjt, l,y) -pjt - s, l,y) ^ ^_ t hjy/jt-s)) y‘^s hjy/jt - s)) 

pit,i,y) ~ t-s hiy/t) 6tit-s) hiy/t) 


Observe that ^ < 1/3 whenever s < t/A. Thus, using (|2.ip and < 1 + 62/6 for 

2 G (0,1/3), we obtain 


t hjy/jt-s)) 
t - s hjy/t) 


- 1 < 


< 


ih (“p iw^) 

6ys s 

5(t — s)2 t — At'^/y'^ 



+ 


s 


t-s 



/ 32 t \ 
Vl5y y2 - At) ' 


Moreover, the monotonicity of /o(' 2 ) gives 

Pit,l,y)-pjt-s,l,y) ^ y^ / 32 ^ ^ 1 

pitA,y) ~ t‘^^ V15y y^-A:tJ Gtjt-s) 

>yl ^ 

“ 15y y^-4:t) 156^2’ 


since y > 52 and y^ — At > 12t. Hence we have showed that 

Pit,l,y) -pit-s,l,y) _ y^ 

pitA,y) ^ ^ 
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whenever s < 4:t‘^jy‘^. Thus, taking into consideration the estimates of Q'xiC'S) for s < 2x (see 
m) , we can write 


V f 1 
{x - 1)2 


exp 


(x-lf 


2s 


2 „,2 poo 


ds = 


(x - 1 )^ y 

y/2 t2 


u 


V2 


r -1/2, 


(x — 1)" 

4x 


where r(a,z) denotes the incomplete Gamma function. Since T{a,z) ^ as z ^ 0'^ and 
r(a, z) ~ 2““^ exp {—z) if 2; —>■ 00 we obtain the following 

,2 




r 

f2 


(^(x - 1) A v^exp ^ . 


On the other side, the estimates of for s > 2x (12.81) gives 


To ~ In X • 


y2 r4t2/y2 

J 2 X 


1 


2x ln(s)ln(s/x) 


exp 


2 s 


ds 


(3.3) 


(3.4) 


and obviously 


hfi/y^ lnsln(s/x; 


■ exp 


("-TO, > 21 ^ 


2 s 


exp 


IP 


=-1/4 


> 


y2 ln{4:t‘^/y'^)ln{4:t‘^/{xy‘^)) 2 y"^ In^{3t/y) 

In a similar way, dropping the exponential term in the integral in (13.4p and estimating the 
logarithms by its values at 2t^ jy^ lead to 


/ 

J2i 


AP I y^ 


ds 


< 2 


I 


ds 4f2 

< 


hp/y^ ln(s)ln(s/x) J 2 P/y^ In^s y^ ln‘^{2t‘^/y"^) y^ln{3t/y)' 

where the hrst inequality follows from s^ jx^ > 2t‘^ j(x‘^y‘^)s > s. Collecting all together we have 
just obtained that 


Inx 


In {3t/y) 


(3.5) 


(y - 1 )' 


2t t — s 


(3.6) 


Finally, to deal with T 3 , observe that m gives 

J ^ p(t,l,y)-p(t-s,l.y) > 1 _ 

p{t,l,y) t-s 

If s G jy‘^,tj2) then t/{t — s) <2 and — It ^ (since y > 52). Consequently, 

the right-hand side of (|3.6I) is bounded from below by 1 — 2e“^ > 0. On the other hand, since 
> 2{t-s) s > f /2 and the function 

{y - 1 )^ t 

At t — s 

is decreasing on {t/2,t), we get the lower-bounds of the form 


s i-A 


exp 


pitA,y) - pjt - s,i,y) 
p{t,i,y) 


t — s 


>1 — 2 exp — 


2t 


>l- 2 e"^> 0 , se{t/2,t). 


This together with the estimates of (?a:,i(s) for large s give 


/ qx{s)ds PS Inx / -- 


1 


exp 


(x — 1)2 

2^ 


ds. 


14 ^ 2 iy 2 s ln(s -|- x) ln(l -|- s/x) 

and since 2 s>s-|-x>s, 1 -|- s/x > s/x > ^/s and (x — l)^/( 2 s) < (x — l)^f^/( 8 y^) < 1/8 on 
{At^/y‘^,t) we have 

rt 


A 2 « Inx 


L 


ds 


4 p/y 2 sln^ s 


= Inx 


1 


1 


In( 4 t 2 /y 2 ^ ln( 2 t) 


(3.7) 
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Now, by compare (13.3p and (13.Sp . we can observe that A 2 dominates Ai. 
Hence combining together (|2.9I) . (13.5p and (13.7p we obtain 


1 


Inx 


/■oo fill 

qx,i{s)ds + A 1 +A 2 + A 3 , Inx + ln^{3t/y) ln(3t/y) ln(2t)y '' ln(3t/y) ’ 

where the last estimates holds since xy/t < 1 and consequently l/ln(3t/y) dominates 
l/ln^(3t/y). Finally, since (see (13.21) 1 

,2 


p{t,x,y)-p{t,l,y) . , .xr 

this part is also dominated by lnx/ln(3t/y) and the proof is complete. 


□ 


4. Estimates for xy/t large 

In this Section we focus on proving (|1.4I) . The main ideas of the proofs in this case come 
from [T], but for convenience of the Reader we adapt the arguments to the case p = 0 and 
present them in condensed form in the following two proposition related to the lower and upper 
bounds respectively. Similarly as in the previous section, we assume that y > x > I, but now 
we consider xy > t . 


Proposition 3. There exist constant Ci > 0 such that for every x,y > I and t > 0 we have 

{ x - l ){ y - l )\ 1 f { x-yf 


Piit,x,y) > Cl 1 A 


t 




exp 


2 t 


(4.1) 


whenever xy > t. 


Proof. For any Borel set A C (l,oo), using the absolute continuity property (12.5D and simply 
estimating the exponential term by 1, we obtain 


Pi { t , x , y ) m { dy ) = 


Ai/2) 




> t,Rt e A;{ — ] exp Q j 

> (Ri)-C2] 

= ■ p^^^‘^\t,x,y)rndl‘^\dy). 

Ja 

Since m ^^^‘^\ dy ) = y ■ m { dy ) it follows that 

Pi { t , x , y ) >(x yy /^ p [^/^\ t , x , y ). 


*d^ 

R? 


(4.2) 


Thus, we conclude by using the explicit formula for p^^^‘^\t,x,y), which is just the Dirichlet 
kernel for the classical Laplacian multiplied by the factor {xy)~^ (see for example [2] for more 
details) 


V^xy 


exp 


2 t 


exp 


, (x-l)(y-l)\ 1 1 

1 A - — - - -p exp 

t J xy^/t 


2 t 

jx - yf 

2 t 


□ 


Proposition 4. There exist constant C 2 > 0 such that for every x,y > I and t > 0 we have 


Pi{t,x,y) < C2 



l)(l/-l) ^ 


1 


y/xyi 


exp 


{x-yf \ 

2 t ) ’ 


whenever xy > t. 
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Proof. Notice that if the space variables are bounded away from the boundary, i.e. y > x > 2, 
then we have 

{x - l){y - 1) ^ ^ 1 

t - 4t 4 

and consequently the result follows from pi(t,x,y) < p{t,x,y) and the asymptotic expansion of 
Io{z) for large arguments (see (12.2^ 1. which implies the estimates of p{t, x, y) of the desired form 
for xy > t. Moreover, the result for small t, i.e. t < 4, follows in a similar way as in Proposition 
[3l Indeed, since 

11 

on using the absolute continuity property (12.51) . we obtain 

Piit,x,y) < {xy)^/^p[^^‘^\t,x,y) 


and the estimates of p^^‘^\t,x,y) gives the result. 

Thus, it is enough to consider 1 < a: < 2 and t > 4. Notice that since xy > t, we have 
y > A/x > 2. Moreover, for 0 < s < t, the following inequalities hold 


/< 


y 


J — s 

and consequently, we can write 


>/o(- 


> 


t — S — s’ 

loiy/t) /■* 1 f y^ + 1 


ri{t, x,y) = qx,i{s)p{t - s, 1, y)ds > (~ 2^t - s) ) 


Now using the following lower bounds for qx^i (see Lemma 4 in [3]) 

(x - 1)2 


X — 1 1 

gx.i(s) > ^=^exp 
V27rx 


2s 


we arrive at 


ri{t,x,y) > 


f 


x-i io{y/t) 

\/2'KX 'Jl Jo y't — S 5^/2 


exp 


(x — 1)^ 
2^ 


s > 0, 


exp 


Substituting ri; = l/s — 1/tin the integral above, we can reduce it to 


,- 1/2 


exp 


(x - 1)2 X 


2t 


w 


- 1/2 


exp —w 


(x — 1)" 


exp 


y^ + i 

2 {t - s) 

y^ + ii 

2/2 yj 


ds. 


dw. 


One can recognize that the above-given integral appears in the integral representation of the 
modified Bessel function Ki/ 2 iz) = j{2z)e~^ and consequently, it can be computed explicitly 
(see for example formula 3.471.15 in [9]). Finally, we get the estimates of the ri{t,x,y) in the 
following form 

(x - 1 + \/pTT)2\ 


u io{y/t) 

ri{t,x,y) > ——-exp 


'X t 

On the other hand, because of (|2.1I) . we have 

pit, X, y) < exp 

and consequently, from (14.3p and (14.4p . we get 

(x - yf 


jx - yf 
2 t 


2 t 


exp ( 


(4.3) 


(4.4) 


Piit,x,y) < 


2 t 


, y 

exp -- 


1 - 


exp 


(x- l)(vV + l + 2/- 1) 


2 t 


Now observe that 

1 f ix-l)iy/y‘^ + l + y-l)\ x-1 f (x- l){y/y‘^ + l + y - 1) 

^ ^ “P (- 2t -j - ^ “P (- 2t - 
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Estimating the exponential term above by 1 and taking into account that for 1 < x < 2 < y 
and xy >t we have 

x-1 ^ {x-l){y-l) t ^ _ {x - l){y - 1) 

y/x t y/x ■ {y - 1) ~ y/x t 

together with the obvious estimate 


1-;= exp 

y/x 


(x- l)(x/y^ + l + y- 1) 


2 t 


< 1 , 


we arrive at 


p,(t, P, p) < _ (^1 A- - -j exp 


jx - yf 
2 t 


, y 

exp -- 


Finally, the asymptotic description of loiz) at infinity (|2.2D ensures that there exists constant 
Cl > 0 such that Io{y/t) < ciy/t/y ■ and consequently 

{x-yf' 




2 t 


where we have additionally used the fact that 1 < lo^xy/t). Hence the proof is complete. 


□ 
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